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ANOMALOUS NODAL COUNT AND SINGULARITIES IN THE 
DISPERSION RELATION OF HONEYCOMB GRAPHS 

RAM BAND, GREGORY BERKOLAIKO, AND TRAGY WEYAND 


Abstract. We study the nodal count of the so-called bi-dendral graphs and show that it 
exhibits an anomaly: the nodal surplus is never equal to 0 or /3, the first Betti number of the 
graph. According to the nodal-magnetic theorem, this means that bands of the magnetic 
spectrum (dispersion relation) of such graphs do not have maxima or minima at the “usual” 
symmetry points of the fundamental domain of the reciprocal space of magnetic parameters. 

In search of the missing extrema we prove a necessary condition for a smooth critical 
point to happen inside the reciprocal fundamental domain. Using this condition, we identify 
the extrema as the singularities in the dispersion relation of the maximal abelian cover of 
the graph (the honeycomb graph being an important example). 

In particular, our results show that the anomalous nodal count is an indication of the 
presence of conical points in the dispersion relation of the maximal universal cover. We also 
discover that the conical points are present in the dispersion relation of graphs with much 
less symmetry than was required in previous investigations. 


1. Introduction 

Quantum graphs and discrete graphs play a signihcant role in numerous recent investi¬ 
gations in mathematical physics. On one side, their importance is apparent in the helds of 
quantum chaos and spectral theory |KS991 IGS061 IKucOSl IBK13] . On the other side, they are 
applied to model various quasi one-dimensional physical systems such as photonic networks, 
nanostructures, and waveguides |Kuc021 IBK131 IPosl2j . The current work exposes a curious 
link between what has so far been considered a more theoretical aspect of spectral theory 
on graphs, the nodal count of eigenfunctions, and a phenomenon with wide-ranging applied 
consequences, the existence of Dirac points in the dispersion relation of periodic structures. 

The study of zeros of a quantum graph’s eigenfunctions is a fertile area of research with 
many questions which are still open. Some recent results include bounds on the number of 
nodal domains |PPA096l. IGSWOdl ISchObf IBerOS] , specihc formulae for some classes of graphs 
|BBS12] . variational characterizations |BBRS1^ IBW13] . and inverse problems |Banl4j . To 
highlight just one key result, it was proven in a sequence of works by different authors (see 
|Ber081 IBBRS12] and references within) that the number of zeros of the n-th eigenfunction 
obeys the bounds 

0 < 0n - (!!■ - 1) < /5, 

where {3 is the hrst Betti number of the graph (intuitively, the number of cycles). A natural 
next question is to consider the distribution of the values — {n — 1)}^^ One would 
expect, for example from the magnetic variational characterization of the “nodal surplus” 
— (r — 1) |Berl31 IBW13j . that all the integers between 0 and [3 appear inhnitely often 
as the nodal surplus of any graph with Betti number f3. However, in the current work we 
show that there is a family of graphs for which this distribution is not supported on the lower 
and upper bounds (0 and (3). Those are the graphs which are obtained as two copies of a 
tree graph glued together at their corresponding leaves. Hence they are called bi-dendral 
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graphs (see figure 3.1[ a) for an example). When the underlying tree is actually a star graph 
(i.e., a graph with one central vertex connected to all other vertices which are degree one), 
we call the resulting bi-dendral graph a mandarin graph. We prove that the above nodal 
count anomaly implies the presence of special singularities in the dispersion relation of the 
abelian cover version of the mandarin graph. This abelian cover is a periodic infinite graph, 
also known as the honeycomb graph; it is a tiling of the plane by congruent hexagons whose 
parallel edges are of equal length (figure me)). 

Periodic infinite quantum graphs have been fruitfully used to model diverse physical sys¬ 
tems such as photonic crystals |KK02j . graphene |KP07j . and its allotropes |DK13j . The 
Floquet-Bloch theory (see chapter 4 in |BK13j i reduces the problem of determining the 
continuous spectrum of a periodic graph to the study of a parameter-dependent opera¬ 
tor on a compact graph. The parameters can be interpreted as magnetic fluxes through 
the graph; the (now discrete) spectrum as a function of these parameters is called the 
dispersion relation. The points where two sheets of the dispersion relation touch are of 
particular interest, as many physical properties of the material are related to the loca¬ 
tion of these points (so-called Dirac points) and the structure of the bands in their vicinity 
|CNGP'*~09l INovlll IFW12[ IKatl2( IFW14j . The current work characterizes the location of the 
Dirac points and the structure of the corresponding eigenstates for periodic graphs, which 
are the abelian covers of the mandarin graphs. This is particularly important as one member 
of this family (the 3-mandarin) is exactly the well-studied hexagonal lattice which models 
graphene, although with much reduced symmetry. The presence of the Dirac points is the ex¬ 
planation for the anomaly in the nodal count (see section]^ below); conversely, the anomalous 
count is an indication of existence of singularities in the dispersion relation. The connection 
between the two is provided by the “magnetic-nodal” theorem |Berl31 ICdV131 IBW13] . 

The next section introduces Schrodinger operators on quantum graphs and section 
presents the main results of the paper. Section contains the proof of the anomaly of 
the bi-dendral graph nodal count, to which we arrive by establishing new eigenvalue inter¬ 
lacing results. Section studies the dispersion relation of general periodic graphs and then 
discusses the abelian cover of the mandarin graph and makes the connection with its nodal 
count. Finally, in appendix]^ we adapt the results of sectionto discrete graphs. 


2. Schrodinger operators on quantum graphs 

We start by defining a quantum graph, following the notational conventions of jBK13] . 
which also contains the proofs of the background results used in this section. 

Let F be a compact metric graph with vertex set V and edge set E. Let iL^(F, C) be the 
space of all comp lex-valued functions that are in the Sobolev space H^{e) for each edge, or 
in other words 


//'■(r.C) = 0ff'-'(e). 

eeE 

Consider the Schrodinger operator with electric potential g : F —)■ M defined by 


fl 

dx^ 


+ g/, 
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acting on the functions from C) satisfying the 5-type matching conditions 

{ f\x) is continuous at n, 

^ = Xvf{v), X,; e M 

ee-B„ ® 

at all vertices v & V. Here the potential q{x) is assumed to be piecewise continuous. The 
set is the set of edges joined at the vertex v] by convention, each derivative at a vertex is 
taken into the corresponding edge. We denote by Xe the local coordinate on the edge e. On 
vertices of degree more than one we usually take Xv = the so-called Neumann condition. 
Non-zero Xv will be used on vertices of degree one, where we also allow the Dirichlet condition 
f{v) = 0, which is formally equivalent to Xv = oo. We note that a Neumann vertex of degree 
two can be absorbed into its neighboring edges unifying them both to a single edge (whose 
length equals the sum of both), without changing the graph spectral properties. At vertices of 
degree two we will sometimes be using the so-called anti-Neumann vertex condition. Namely, 
for a vertex v of degree two, which is connected to the edges e+ and e_, the anti-Neumann 
condition is 


(2-2) /|,_ (u) = - (u), /'L_ (u) = (u), 

where the direction of the derivative is taken from the vertex v into the edge e±. Note that 
the anti-Neumann condition does not fall in the class of vertex conditions (2.1). 

The operator is self-adjoint, bounded from below, and has a discrete set of eigenvalues 
that can be ordered as 

-5i < A 2 < ... < A„ < ... . 


Throughout the paper we will frequently assume that the graph’s eigenvalues and eigenfunc¬ 
tions are generic in the following sense: 


Definition 2.1. 

(1) The eigenpair (A, /) is called generic if the eigenvalue A is simple and the correspond¬ 
ing eigenfunction f is different than zero on every vertex. 

(2) Both A and f in a generic pair are also called generic. 

(3) A quantum graph is generic if all of its eigenpairs are generic. 


The conditions which guarantee the graph’s genericity are discussed in |Fri05( IBL14j . We 
will count zeros only for generic eigenfunctions. 


Definition 2.2. 

(1) Let fn be a generic eigenfunction corresponding to the n-th eigenvalue, An (counting 
with multiplicities). Denote by the number of its internal zeros. Namely, we do 
not include the Dirichlet vertices, if they exist, in the count. 

(2) The quantity an ■= (fn — {n — 1) is called the nodal surplus. 


It was recently discovered that the graph’s nodal count is closely related to properties of 


the magnetic Schrodinger operator on the graph. This connection is described in theorem |2. 3 
below. The magnetic Schrodinger operator on T is given by 




f + qf, feH\T,C), 
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where the magnetic potential, A{x), is a one-form (namely, the sign of A{x) changes with the 
orientation of the edge). The 5-type boundary conditions are now modihed to the following 
at all vertices v eV\ 


{ /(x) is continuous at u, 

^ (^(^) “ = Xvf{v), Xv e M. 

Let /3 = |£'| — |l/|-|-lbe the hrst Betti number of the graph T, i.e. the rank of the graph’s 
fundamental group. Informally speaking, f3 is the number of “independent” cycles on the 
graph and hence is zero if the graph is a tree. Up to a change of gauge, a magnetic held on 
a graph is fully specihed by /3 huxes ai, 02 ,..., dehned as 


aj = (p A{x)dx mod 27r, 

where {tj} is a set of generators of the fundamental group. In other words, magnetic 
Schrodinger operators with different magnetic potentials A{x), but the same huxes (ai,..., ajs), 
are unitarily equivalent. Therefore, the eigenvalues \n{H^) can be viewed as functions of 
CK = («!,..., ajs). The connection between this function and the nodal count is explicated 
in the following theorem. 

Theorem 2.3. Let {Xn, fn) be a generic eigenpair of with A corresponding to some flux 
Q:*G{0,7r}^. Then ct* is a non-degenerate critical point of the function A„(q:) and its Morse 
index is egual to the nodal surplus, an- 

Remark 2.4. The above theorem is proved in |BW13] for cx. = (0, ...,0), extending to 
quantum graphs an earlier result on discrete graphs |Berl3j (also see |CdV13] ). One may 
repeat the same proof from |BW13j for all other a G {0, tt}^, noticing the following: 

(1) Introducing a vertex of degree two with an anti-Neumann vertex condition on some 
edge results in an operator which is unitarily equivalent to the original operator with 
a magnetic potential on this edge integrating to vr. 

(2) The anti-Neumann condition keeps the operator self-adjoint and does not change 
either of the quadratic forms used in section 4 of |BW13j : the only change needed is 
the additional condition /|^ {v) = —f\^ (v) imposed on the domain of the quadratic 
form. 

(3) The eigenfunction of a tree graph has n — 1 internal zeros and this holds even if 
some anti-Neumann conditions are imposed on the graph (relevant for the proof of 
theorem 3.3, part 1 in |BW13] L 

As a corollary, we get an earlier result (see [BernSl IBBRSl^ and references therein), 

(2.3) 0<an<f. 

The theorem above makes one wonder whether the whole range of integers from 0 to /? is 
covered by cr„ for any graph. We will show that there is a family of quantum graphs for 
which the answer to this question is negative. These graphs, to be discussed next, and their 
spectral properties are the focus of this paper. 
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Figure 3.1. An example of a bi-dendral graph (a), the symmetry axis is 
indicated by the dashed line; a 3-mandarin graph (b) and its maximal abelian 
cover (c). The abelian cover graph has honeycomb structure although less 
symmetry than a regular honeycomb lattice: only the parallel edges have the 
same length. 


3. Main results 


3.1. Anomalous nodal count. Consider two copies of a tree graph with d leaves (vertices 
of degree one) which are glued together by identifying corresponding leaves of both. We 
impose Neumann conditions at all vertices and call the resulting quantum graph a bi-dendral 
graph. When the underlying tree graph is chosen to be a star graph, we obtain a graph 
consisting of two vertices and d edges connecting them (recalling that Neumann vertices of 
degree two can be absorbed into an edge). We will call such a graph a d-mandarin graph 
(see hgure 3.1). A bi-dendral graph has an obvious symmetry axis which passes through 


all points arising from gluing the leaves. We call the edges which cross this symmetry axis 
the middle edges. We extend the dehnition of a bi-dendral graph by allowing for vertices of 
degree two with anti-Neumann condition to be present on those middle edges. Finally, we 
will be assuming the graph edge lengths are chosen such that either a particular eigenvalue 
is generic or the graph is generic (see dehnition 2.1). 


Theorem 3.1. Let T he a bi-dendral graph with d middle edges, such that a > 0 of them have 
anti-Neumann conditions imposed at an intermediate point. If the eigenpair is generic 
and n > 1 or a > 0, then the nodal surplus an satisfies 

(3.1) l<an<(3-l = d-2. 

Remark 3.2. Compare (3.1) with the general bounds of the nodal surplus, (2.3). 


Consider the eigenvalue An(Q:) as a function of ck = (oi,... , 0 /?), the total huxes through 
the cycles of the graph (for some pre-hxed basis of the fundamental group). It is well-known 
that the points cx G {0, vr}^ are critical points of the function Xn{ct). they are the hxed points 
of the symmetry transformation cx i—)■ —cx (mod 27r). Henceforth we will refer to these points 
as symmetry points. 

Theorem |2.3| states that the nodal surplus of the eigenfunction at a symmetry point is equal 
to the Morse index of A„(q:) at this point. We now infer from theorem 3.1 that none of the 
symmetry points are extrema of the eigenvalues An,(Q:) (with the exception of the eigenvalue 
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Figure 3.2. First four bands of a 3-mandarin graph, as functions of /3 = 2 
magnetic fluxes. The vertical axis is k = a/A; in this scale the bands have 
comparable sizes |BB13j . This choice of scale also explains the conical point 
at the bottom of the lowest band (in contrast, Ai(q:) is differentiable at 0). 


Ai(0,..., 0), which is always a minimum). Since the A„(q:) are continuous functions on the 
compact torus (—vr, tt]^, the extrema must be achieved somewhere else. 


3.2. A condition for an internal critical point. Searching for the “missing” extremal 
points of the functions A„(q:), we need to look in the bulk of the fundamental domain of the ck- 
space. Critical (extremal) points that are not one of the previously identihed symmetry points 
will be called internal critical (extremal) points. A common (and mistaken) assumption is 
that there are no internal extremal points; it was pointed out in |HKSW0~7] that such points 
may occur and explicit examples were constructed in |HKSW07l lEKWlO] . Moreover, in 
the case of mandarin graphs their occurrence is simply unavoidable; the extrema must be 
achieved somewhere. 

Our search for extremal points is aided by proving a necessary condition that simple 
internal critical points must obey. 


Theorem 3.3. Let T be a graph with first Betti number equal to fi and denote by ex & (~'^5 
the total fluxes through some choice of cycles of the graph that form a basis of its fundamental 
group. If the eigenvalue A„(q:) has an internal critical point ex* (i.e., cx* ^ {0,7r}^j and 
A„(q:*) is simple, then the eigenfunction corresponding to A„(q:*) is equal to zero at some 
vertex of the graph. 


3.3. The missing extrema are found at touching bands. The necessary condition given 
in theorem 3^ allows us to show that an internal critical point of a mandarin graph must be 
degenerate. Calling the graph of the function An(cK) a band, we get the following: 


Theorem 3.4. For a generic mandarin graph, all extrema of Xn{(x), apart from the minimum 
of Xflex), are achieved at points where two bands touch. 
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Figure |3.2 illustrates the result of the theorem: each pair of consecutive bands touch at 
two symmetrically located points. 

It is well-known that the magnetic spectrum of a graph as a function of the fluxes coincides, 
by Floquet-Bloch reduction, with the dispersion relation of the abelian cover of the graph. 
Here the phases a take the meaning of the quasi-momenta of the Bloch functions. For a 3- 
mandarin graph there are two a-parameters, and a degeneracy in two parameters is typically 
a conical point |Arn781 app. 10]. The maximal abelian cover of a 3-mandarin is a honeycomb- 
like graph, where only the parallel edges are required to have equal length (see hgure 3.1). 
The above chain of results shows that typically there are conical (“Dirac”) points in the 
dispersion relation of such a graph between each pair of adjacent bands. 


4. Interlacing and nodal count for mandarin graphs 


4.1. Eigenvalue interlacing. We begin working towards the proof of theorem 3.1 by es¬ 


tablishing some eigenvalue interlacing results, which are heavily based on the inequalities of 
|BK13i sect. 3.1]. The following lemma was proved as part of lemma 4.3 of |BBS12] . 

Lemma 4.1. Let T be a compact connected graph and let A = An(F) he a simple eigenvalue 
of T with the corresponding eigenfunction f. Let Fc be a graph which is the union of k + 1 
connected components such that the following are true: 

(1) F can be obtained from Fc by k operations of gluing a pair of vertices together and 
adding the parameters x of their 6-type conditions. 

(2) The function f restricted to any of the components of Fc is an eigenfunction of that 
component. 


Then 

(4.1) 

(4.2) 


An—l(Fc) < A < An+fc+l (Fc), 

An(rc) = An+l(Fc) = . . . = An+fe(Fc) = A. 


The next lemma we formulate in somewhat greater generality than that which will be 
necessary in subsequent derivations. 

Lemma 4.2. Let T he a graph with a Neumann vertex v of degree d whose removal separates 
the graph into d disjoint subgraphs. We denote its edge set by E^. Let r < d be a non-negative 
integer. For a subset Ed of E^, with \Ed\ = r, define F^^ to be the modification of the graph 
F obtained by imposing the Dirichlet condition at v for edges from Ed and leaving the edges 


from Ey \ Ed connected at v with the Neumann condition (see figure f.l for an example). 
Then 


(4.3) 


:i..-i(r) < mill \(rEo) < A„(r) < max A„(rEo) < A„+i(r). 

\Eu\=r \EB\=r 


Proof. Without loss of generality we will assume that the eigenvalue An(F) is simple. Indeed, 
one can resolve multiplicity by an arbitrarily small perturbation of the edge lengths |Frin5] . 
Since the n-th eigenvalue is a continuous function of the lengths |BK13l thm 3.1.11], if the 
(non-strict) inequality is true for simple eigenvalues, it remains true when passing to the limit 
of zero perturbation. 

The outside inegualities: Start from the graph F and disconnect any r edges keeping them 
joined together with the Neumann condition and denoting the resulting graph by F'. From 
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(a) (b) 

Figure 4.1. The original graph T, (a), and the resulting graph T^;^, (b), 
when the set Ed comprises of the r = 2 edges joining v from the left (shown in 
thicker lines). Vertices with Neumann conditions are shown as empty circles, 
and Dirichlet conditions are indicated by hlled circles. 


jBK13I thm 3.1.11] we get 

An-i (r) < (F) < A„+1 (F) < A„+1 (T). 

Now, changing the vertex condition from Neumann to Dirichlet disconnects the r edges and 
results in A„ (F) < An (Fe^) < An+i (F') |BK131 thm 3.1.8]. 

The lower bound for the maximum: We describe the choice of Ed that will fulhll the 
inequality. We start with the the n-th eigenfunction / of F (normalizing so that f{v) > 0). 
We now choose the r edges e with the largest values of fe{v) to be in Ed and get 

E /'w<o. 

eeE\Eo 

Consider / on the graph F'^^ obtained by disconnecting the edges from the set Ed at the 
vertex v. We supply F^^ with conditions that ensure that / is an eigenfunction. Namely, 
the function / satishes the 5-type condition at the vertex v with some coefficient Xv F 0. At 
the new vertices Ui, V 2 etc. (see hgure 4.1) some 5-type conditions are also satished. The 
sum of the parameters of those 5-type conditions together with Xv is zero. From lemma 4.1 
we have 

A„(riJ = A„(r). 

Now we modify the conditions of F^^, increasing Xv to zero and the 5-type parameters at the 
new vertices of degree one to oo (Dirichlet condition), obtaining the graph Feo- All these 
operations increase the eigenvalues |BK13l thm 3.1.8]; in particular, 

A„(r£„) > A„(r'jj = A„(r). 

The upper bound for the minimum: In this case, we choose the r edges e with the smallest 
values of /e(n) to be in Ed and get 

E /;(«)> 0, 

e&E\Eo 

and, therefore, Xv F 0. From lemma M we have 

A„+.(rij = A„(r). 
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Now we decrease Xv to 0, resulting in 

^n+r(^Eo) — -^ri,+r(r£;3)5 

and then increase the parameters of the new vertices of degree one to oo, obtaining 

An(rE^) < Xn+r{^Eo)- 

The desired result is an obvious chaining of the obtained inequalities. □ 


Remark 4.3. The condition in lemma |4^ that removing a vertex of degree d separates the 
graph into d disjoint subgraphs is necessary for the lower bound in (4.3). This is used in 
the proof above when applying lemma 4T, and indeed there exist examples where the bound 
does not hold when the condition is not satished. 


Remark 4.4. In the case of a graph which is “spherically-symmetric” around vertex v (i.e. 
the d disjoint subgraphs obtained by removing the vertex v are identical), we immediately 
get 

An(r_E£,) = '^n(r), 

for all n and Eq . A special case of this is proved and used in a recent work of Demirel-Frank 
[DTT5] . 

Lemma 4.5. Let Ten be a tree graph with zero potential, Dirichlet conditions at t of its 
leaves, and Neumann conditions at all of its other leaves and at all of its internal vertices. 
Let Tnd be the same graph with Dirichlet conditions changed to Neumann and vice versa at 
all leaves. If k + t — 1 > 0, then 

(4.4) Xk{TEN) < Afc+t-l(rAri5). 

Remark 4.6. Naively applying the standard interlacing results (see, e.g. |BK131 thm 3.1.8]) 
while changing the condition at each leaf, we would get the weaker result 

Xk{T en) < Xk+t{TND)- 

Remark 4.7. Taking the inequality above with t = 0 gives an inequality similar to that 
known for domains in M'^. In general, if hi C is a domain whose Neumann spec¬ 
trum is discrete, then A^^^(r2) < A[,^^]^(r2) where the superscript {N)\{D) denotes the 
Neumann\Dirichlet spectrum |Fri91( iFilOdj . 


Proof. We start by introducing the spectral counting function 

7 V(F;A) :=|{n|A„ (r)<A}|, 

which we will use to rephrase eigenvalue inequalities. In particular, we have the following 
equivalence: 


(4.5) Afc (Fi) < Xk+s (ha) for all fc N (F 2 ; X) < N (Fp A) + s. 

To go from left to right we observe that if A^(Fi; A) = n, then X < A^+i(Fx) < A„_|_5 ^i(F 2 ), 
by dehnition and left inequality, correspondingly. Applying the dehnition of N once again, 
we get A^ (F 2 ; A) < n + s. The left follows from the right by the substitution X Xk+s ( 1 ^ 2 )- 
Using the equivalence above, we write (4.4) as 

(4.6) N {T^e-, X)<N (F^^,; A) + f - 1, 


which is what we now prove using induction on the number of internal vertices of Ten- The 
starting point is to notice that the lemma holds for intervals (no internal vertices) with either 
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Dirichlet or Neumann conditions at each of their endpoints. Denoting intervals with different 
types of boundary conditions by IdD) Inn, and Idn = Ind, we have 

/TT \ 2 

Afc {Idd) = Afc+1 {Inn) = 

Afe {Inn) = Afc-i {Idd) = (y “ ^)) 

Afc {Idn) = Afc {Ind) = ( 2 fc — 1 ) j , 

where I is the interval length. Therefore the lemma holds for an interval, and for such graphs 


we even get an equality in (4.4). Assume that the lemma holds for all tree graphs which 


have no more than M internal vertices. Let r£)jv be a tree graph which satishes the lemma’s 
conditions and has M + 1 internal vertices. Let v be an internal vertex of L^jw with degree d. 
Note that v has a Neumann condition as it is an internal vertex. We start with Xk (L dn) and 
apply lemma 4.2 with r = d — 2 on the vertex v (choosing Ed which results in the maximal 


eigenvalue). Denoting the resulting graph by L', we have from lemma 4.2 that 
(4.7) Afc (rzuv) < Afc (D), 


or equivalently from (4.5), 

(4.8) iV(r';A) <iV(rB^;A). 

The resulting graph, T', is a union of d — 1 tree graphs, T' = Note that d — 2 of 

these graphs have a Dirichlet leaf vertex which originates from v, and therefore these graphs 
have less internal vertices than Tdn- In one graph out of the T^dn^ however, v becomes a 
Neumann vertex of degree 2 and can be absorbed into the edge. Therefore this graph also 
has strictly less internal vertices than Tdn- Denote by ti the number of Dirichlet leaves of 


T^dn 


d-l 

E 

2=1 


ti — f T d — 2, 


since we added d — 2 Dirichlet vertices by splitting the vertex v. 
We get by the induction assumption 


4.9) 


iv(D>D;A) <iv(rgyA) +«,-! 


where T^^^ is the same as T^^jy, but with Dirichlet and Neumann conditions at all leaves 
interchanged. 

Denoting T" = U)ti and using the additivity of the spectral counting function we get 


d-l 


d-l 


(4.10) N (T"; X) = J2^ [^nd-^ ^ E X) + U - I)= N (D; A) + f - 1. 

2=1 2=1 

All that remains is to note that Tat^) is obtained from T" by gluing all vertices which origi¬ 
nated from V and are equipped (in T") with Neumann conditions; such gluing increases the 
eigenvalue [BK131 thm 3.1.11], so A^+i-i (T") < Xk+t-i {Tnd) and therefore 

(4.11) iV(r^z);A)<iV(r";A). 

Combining inequalities (4.8),(4.10),(|4.11[), we get (4.6), which proves the induction step. □ 
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Figure 4.2. A bi-dendral graph with an anti-Neumann condition on the 
fourth middle edge, indicated by a cross. The symmetric eigenfunctions are 
eigenfunctions of the top tree graph on the right; the antisymmetric ones are 
the eigenfunctions of the bottom tree graph. Neumann conditions are shown 
as empty circles and Dirichlet conditions as hlled circles. 


4.2. Application to nodal connt. The interlacing results of the previous section allow 
one to prove that bi-dendral graphs have anomalous nodal counts: the nodal surplus of a 
bi-dendral graph with d middle edges is never /3 = d — 1 or 0 (with the exception of ai which 
is 0 for all graphs with Neumann conditions). We will prove this fact for a slightly more 
general setup. 

We consider the bi-dendral graph, some of whose middle edges can have an intermediate 
point V at which the anti-Neumann conditions (2.2) are enforced. We can still ask the 
questions about the number of zeros of the eigenfunction. Note that there is usually a sign 
change at the above v but not a zero. It can also happen that f\^_ (n) = (v) = 0, but 

then, on the contrary, / does not change sign at v. The n-th nodal surplus, as before, is 
dehned as the number of zeros of the n-th eigenfunction minus n — 1. 

Placing an anti-Neumann condition at a point on the edge e is unitarily equivalent to 
placing a magnetic potential that integrates to tt on the edge e. In particular, without loss of 
generality we can impose each anti-Neumann condition in the middle of the respective edge. 

We are now ready to prove theorem 3T, which asserts that a generic bi-dendral graph 
satishes 

l<<Tn</3 — l = d — 2, 


n > 1. 


Proof of theorem dA Let a be the number of edges on which anti-Neumann conditions are 
enforced. Since the anti-Neumann conditions are symmetric and can be imposed in the 
middle of edges (their location along the edge does not affect the eigenvalues nor the nodal 
count), the graph still has up-down reflection symmetry (see figure 4.2, left). Therefore, the 
spectrum can be decomposed into two halves, even and odd, consisting of eigenvalues whose 
eigenfunctions are symmetric and those with anti-symmetric eigenfunctions with respect to 
the reflection. 

The symmetric eigenfunctions have Dirichlet conditions at the midpoints of anti-Neumann 
edges and zero derivative at the midpoints of the rest of the middle edges. There is a one- 
to-one correspondence between the symmetric eigenfunctions and the eigenfunctions of the 
upper half of the bi-dendral graph, Teven, with the corresponding boundary conditions (see 


hgure 4.2, right). The antisymmetric eigenfunctions satisfy the opposite conditions and are 


in correspondence with the eigenfunction of the tree graph Todd- Namely, if T is a bi-dendral 
graph, then we have the spectral decomposition Spec (T) = Spec (Teven) U Spec (Todd)- The 
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resulting pair of trees have opposite vertex conditions on the leaves and therefore £t the 


description in lemma 4.5 


Consider the eigenfunction corresponding to Afc(reven), and assume that this eigenvalue is 
generic in the spectrum of F. As Feven is a tree graph, this eigenfunction has k — 1 internal 
zeros (see Remark [^ , part|^ and a zeros at the Dirichlet leaves. Its unfolding to the graph 
F therefore has 2{k — 1) + a zeros. Applying lemma 4.5 twice for comparison with the graph 
Fodd (which has t = d — a Dirichlet leaves), we get 

Afc_(i_|_a_|_i(Fodd) < Afc(F even ) ^ Afc^a_i (Fodd)) 


the second inequality being valid only when k+a — 1 > 0, which is true due to our assumptions 
(/c > 1 or a > 0). The inequalities above are strict as we assumed A^(Feven) to be simple when 
considered as an eigenvalue in the spectrum of F. The position of the eigenvalue Afc(Feven) 
in the spectrum of F is equal to the number of eigenvalues of either Feven or Fodd that are 
smaller than or equal to A*,(Feven)- From the inequalities above, we deduce that this position 
is between 2k — d + a + 1 and 2k — 2 + a. The nodal surplus is then between 

2(fc - 1) + a - (2fc - 2 + a - 1) = 1 

and 

2{k-l) + a-{2k-d + a + l-l)=d-2. 

The calculation for the antisymmetric eigenfunctions is identical modulo the change of a 
to d — a, which has no effect on the final result. □ 


In the two simplest cases of mandarin graphs (we assume no anti-Neumann conditions), 
we know all nodal counts exactly. 


Corollary 4.8. The nodal surplus count of a mandarin with 3 edges and no magnetic fluxes 
is {0, 1, 1, 1,...}. The nodal surplus count of a mandarin with f edges and no magnetic 
fluxes is {0, 1, 2, 1, 2, 1,...}. 


Proof. The nodal surplus of the hrst eigenfunction is always 0 (it does not change sign). The 
nodal count of 3-mandarin follows immediately from theorem 3.1| with d = 3. 

Consider now the 4-mandarin graph. The number of sign changes along every edge of the 
graph must have the same parity: the parity is fully determined by the relative signs of the 
eigenfunction on the two vertices of the graph. In either case the total number of zeros is 
even. From theorem 3.1, the (non-trivial) nodal surplus of the 4-mandarin is either 1 or 2, 
but the even number of zeros forces the alternating pattern. □ 


5. Critical points of the dispersion relation of infinite quantum graphs 

5.1. Dispersion relation and magnetic Laplacian. Let G be an infinite Z*^-periodic 
quantum graph with no magnetic potential. The Schrodinger operator on G acts as 

(5.1) {Hf){x) = - fix) + qix)f{x), 

where qf) G M is a bounded, periodic, piecewise continuous function that represents electric 
potential. The conditions on the vertices of G are assumed to be of (5-type (2.1) and are also 
periodic. 

Choose a fundamental domain W, i.e. a connec ted subgraph of G which has one repre¬ 
sentative of each orbit of acting on G (see hgure [Alfa) for an example). We assume that 


the boundary of W does not include any vertices of G (handling the latter would introduce 
unnecessary notational difficulties). 
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(a) (c) 

Figure 5.1. (a) A Z^-periodic graph with a chosen fundamental domain; (b) 
the chosen fundamental domain with quasi-identihed vertices marked with the 
same shape; (c) the graph D formed by merging the quasi-identihed vertices. 


We call a pair of vertices, 

, ) = d 


Denote by some choice of k generators of Z*^ acting on 

(c^,cj), belonging to the boundary of W, quasi-identified if Sj (cj) = (see hgure 5.1(b)). 
We assume that there is only one such pair for each j = 1,..., fc. Note that this condition 
may depend on the choice of the fundamental domain. 

Let ct = (di, Q! 2 , • • •, cpc) and on the graph W dehne the operator : H‘^{W,C) —)■ 
L2(W, C), which acts as — ^ + q{x) on every edge, along with the conditions 


,5 2 , f(c-) = e“-ncl), 

f'(c-) = -e“V'(c+) 

at the quasi-identihed vertices and the 5-type vertex conditions inherited from the graph G 
at all other vertices of W. 

By Floquet-Bloch theory (see e.g., |Kuc931 IBK13] L we know that the spectrum of G can 
be found by calculating the eigenvalues of and taking the union over all ck in the Brillouin 
zone, the torus (—tt, tt]^ of all possible values of ck. In other words. 


a(H) = U 

q;G(— 7r,7r]^ 

The multi-valued function is called the dispersion relation. 

When the operator H is real, as it is in our case, complex conjugation transforms H°‘ into 
implying that the dispersion relation a{H^) is symmetric with respect to the inversion 
CK I—)■ — CK. The hxed points of this transformation are the vectors a G {0, tt}^ with all entries 
equal to either 0 or tt. We call these vectors the symmetry points of the Brillouin zone. We 
remark that if G has additional symmetries, a hierarchy of symmetry points may appear in 
the Brillouin zone (see |BC14j for an example). However, in this work we reserve the term 
for the vectors from {0, tt}^ only. 
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The relation to a magnetic operator on a graph is explicated by the following construction. 
Beginning with the fundamental domain W, glue the vertices of each quasi-identi hed pair 
(c^,c~) to form a new vertex cj and denote the resulting graph by D (see hgure 5.1'c)). 
When connecting vertices, we do not change the edge lengths. The following easy result can 
be found, for example, in |BK13| thm 2.6.1] or |KSn3| [Ruellj . 

Lemma 5.1. The operator is unitarily equivalent to the operator : if^(D,C) —)■ 
L2(D,C), defined as — — iA{x)'j + q{x) on every edge, with the vertex conditions 


g{x) is continuous at v, 
d 


eeE„ 




9{v) = Xv9{v), Xv e 


where A{x) is a one-form on D that satisfies 

rc+ 

I 3 

A{x) mod 27r 


aj — 


for any path on W between cJ and c^. 

The unitary equivalence is as follows. Choose an arbitrary point, p, on W. If g is an 
eigenfunction of H^, then f := ge~'^^ is an eigenfunction of H°‘, where 

i{x) = [ A{x). 


Remark 5.2. The magnetic flux aj is path independent because there is no other magnetic 
potential on our graph. The integral of A{x) around a cycle in W (i.e., a cycle that has no 
magnetic potential) is zero. Since the sign of A[x) (and afi depends on direction, the integral 
around part of a cycle in one direction is equal to the integral around the rest of the cycle 
traversed in the opposite direction. The same applies to the phase ^(x). 

in particular, shows that the eigenvalues of do not depend on the local 


Lemma 5.1 


changes to the choice of the fundamental domain W. Since |/| = and hence /(x) = 0 
if and only if g[x) = 0 , we can assume without loss of generality that eigenfunctions are 
non-zero at Cj, as well as at c^. 

Now we dehne another operator based on the graph W by specifying different condi¬ 
tions at the quasi-identihed vertices. Let 7 = (71, 72 , • • •, 7 a:) and dehne the operator 
: //^(W, C) —)■ L2(W, C), which acts as —^ + q{x) on every edge, along with the 
Robin vertex conditions 

ficj) = -Xjficj) 

at the quasi-identihed vertices and the same conditions as G at all other vertices. It turns 
out that extremal points of A„(q:) manifest themselves as multiple eigenvalues of H''. 


Theorem 5.3. Suppose the infinite periodic quantum graph G has no magnetic potential and 
the fundamental domain W has only one quasi-identified vertex pair in each direction. Let 
A„(q:) have a critical point cx* that is not at a symmetry point of the Brillouin zone (i.e., 
such that a* 7^ 0,7r), and suppose that the eigenvalue An(Q:*) is simple with corresponding 
eigenfunction f. 

Then 










ANOMALOUS NODAL COUNT AND DIRAC POINTS 


15 


(1) is real for all j, 

(2) A = An(Q:*) is a degenerate eigenvalue of , and 

(3) if, additionally, W is a tree, there exists an internal vertex of W, of degree three or 
higher, such that the eigenfunction f is zero at this vertex. 


To prove theorem 5.3 we first collect some auxiliary useful facts. 


5.2. Critical points of the dispersion relation. 


Lemma 5.4. If f is an eigenfunction of the operator ( |5.1| ), then lm{f'{x)f{x)) is constant 
on each edge of the graph. 

Proof. We start by calculating the Wronskian of the functions / and /: 


W{fJ) = nx)f{x)-f{x)f'{x) 


= f'{x)f{x) - f{x)f'{x) = 2ilm{f{x)f{x)). 

On the other hand, both / and / are solutions to the differential equation —y'\x) + {q{x) — 
\)y{x) = 0 on any edge and, by Abel’s Theorem, their Wronskian is constant. □ 

We note that the value of the Wronskian changes from one edge to another. If all vertex 
conditions are Neumann, the Wronskian defines a flow on the graph. This and other facts 
about the Wronskian on graphs can be found in |BW131 Weyl4| . 

In what follows, we will need to differentiate the function An(Q:). This is allowed since 
A„(q:*) is a simple eigenvalue of , and the function A„(q:) is analytic around cx = 
cx* |Kat76] (see also |BK121 IBK13] for a discussion of this fact for quantum graphs). 

Lemma 5.5. Suppose that the infinite periodic quantum graph G satisfies the conditions of 
theorem 5.3 If cx* is a critical point of Xn{cx) and A„(q:*) is a simple eigenvalue, then the 
eigenfunction f of 11°" corresponding to A„(q:*) satisfies 

(5.3) 


/'(T)/(c+) 6 


for all j = 1,2,..., k. 
Proof. We will show that 

(5.4) 


da 


- = -2Im (/'(c+)/(c+)) , 


which directl y im plies (5.3) since ck* is a critical point. 

By lemma 5.1, and have the same eigenvalues so 


dA(Q:) 


daj 


= ^\{cx* + tScXj] 

(Lb 


t=0 


= fMA'+tB) 


t=0 


where 6cXj = (0,..., daj, ..., 0) and B{x) is any continuous function that satisfies 


B{x) mod 271 = 


6ai k = j 
0 k^j 
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Denote by Cj the edge of D which contains Cj. In particnlar, we choose 5aj = 1 and a fnnction 
Bj{x) that is compactly snpported on edge Cj near the point Cj and satishes 


(5.5) 


Bj{x) mod 2tt = 6kj 


where 6kj is the Kronecker delta fnnction. 

Let Qt be an eigenfnnction of norm one corresponding to A (A* + tBj). We define g '■= go 
and get from the Hellmann-Feynman Theorem |Fey39| that 




t=0 


dt 


g,g 


t=0 


One can calculate that 


dt 


t=o 


t=o 

= 2iBj{x )-—h iBj{x) + 2A*(x)Bj{x). 


Since Bj{x) is supported on the edge ej only, we get 
c1A(q:) 


daj 


= J iBj{x)g{x)g{x)dx + J {^2iBj{x)g'{x)g{x)+ 2A*{x)Bj{x)\g{x)fj da; 

[ 9 {x)g{x)^ +2iBj{x)g\x)g{x)+ 2A*{x)Bj{x)\g{x)\‘^'^ da;, 

using integration by parts (the boundary terms disappear due to the support of Bj{x)). 
Continuing, this gives us 


c1A(q:) 


= / Bj{x) -ig{x)g’{x) + A*{x)g{x)g{x) + ig\x)g{x) + Al*(a;)^(a;)^(a;) 


da; 


= / 2BAx)\m. 


-g\x) + iA*{x)g{x)) g{x) da;. 


where mov ing from the hrst line to the second we use the fact that A*{x) is real. By 
lemma 5.1, we know that since g is an eigenfunction corresponding to A(74*), / = ge~'‘^ is an 
eigenfunction corresponding to A(q:*), and therefore 

Im 


{g{x) - iA*{x)g{x)) g{x) = Im /'(x)/(x) . 


However, by lemma 5.4 the latter value is a constant on the edge e^-, and using (5.5) we get 
c?A(q:) 




= -2Im (^f{x)f{x)^ / Bj{x)dx = -2Im (^/'(x)/(x)) . 


□ 


We are now ready to prove the main theorem of the previous subsection. 
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Proof of theorem 5.3. Let / be an eigenfunction of (with a* not at a symmetry point 


of the Brillouin zone), and let 


7 ,- = 


/'(c;) 
/«) ’ 


j = 1,..., fc. 


According to lemma 
and real. 


5.5 


all 7 *’s are real, and consequently the operator H~^* is self-adjoint 


It is easy to see that / is an eigenfunction of . However, / cannot be made real since it 


satishes (5.2) and we assumed that there exists j with a* ^ 0,7r. This is not a contradiction 


only if the real and imaginary parts of / are both eigenfunctions of if'’' , in which case A 
must be a degenerate eigenvalue of if'*'*. 

Furthermore, if W is a tree, we can apply |BK13( cor 3.1.9], which says that if an eigenvalue 
of a tree is degenerate there exists an internal vertex of degree three or higher at which all 
eigenfunctions from the eigenspace vanish. □ 

Theorem 13.31 now follows. 


Proof of theorem 3^. Let T be a graph with hrst Betti number equal to /3 and denote by 
OL G (—TT, tt]^ the total fluxes through some choice of cycles of the graph that form a basis 
of its fundamental group. Therefore, one may cut the graph at f3 positions to make it a 
tree graph. The obtained tree graph serves as a fundamental domain W ( with respect to 
translations) of an inhnite Z^-periodic quantum graph, G, as in theorem 5.3 The statement 
in theorem 3.3 now follows from theorem 5.313), realizing that the quasi-momenta of G are 
exactly the magnetic fluxes of T. □ 


5.3. Zeros and touching bands for mandarin graphs. In light of theorem 5^ we ob¬ 
serve that a special role is played by eigenfunctions that are zero on at least one vertex. We 
now apply this observation to mandarin graphs. 


Proof of theorem 3^. We will now show that if we have an extremum in the dispersion rela¬ 
tion of the mandarin graph, it is due to touching bands. 

Assume the contrary: an extremum of A,i(q:) occurs at ck = ck* and An(o:*) is simple. The 
minimum of Ai always happens at the point a = (0, 0,...) and we exclude it from further 
considerations. 

First we argue that the eigenfunction / corresponding to A„(q:*) must vanish at a vertex of 


the graph. Indeed, if ck* is not a symmetry point, the claim follows directly from theorem 3.3 


On the other hand, if ck* is a symmetry point and is non-vanishing on the vertices, by 


combining theorems |3.1| and |2.3[ we conclude that ck* is a saddle point, which contradicts it 
being an extremum. 

Now, without loss of generality, assume that / vanishes on the top vertex of the d-mandarin 
graph. In addition, / satishes Neumann conditions there. We can shift the magnetic condition 
to the top vertex, resulting in 

(5.6) fi(v) = f 2 (v) = ... = fd(v) = 0, 

(5.7) e‘-/((n) + ... -f + f^(v) = 0, 

the standard Neumann condition at the bottom vertex, and a non-magnetic operator acting 


on the edges. Ignoring, for a moment, condition (5.7), we get a standard star graph with d 


edges and Dirichlet conditions at the boundary vertices. Such a graph, for generic choice of 
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lengths, has eigenfunctions that do not vanish on the central vertex. Hence, we assume that 
/ is not equal to zero at the bottom vertex. 


Apply the top-down (vertical) reflection F to the function / (see figure 4.2 (left)), followed 
by complex conjugation. It is immediate to check that the new function, Ff, satishes the 
same eigenvalue problem with ck = ck* as the function /. It is, however, a different function: 
/ vanishes on the top vertex and does not vanish on the bottom; the function Ff does the 
opposite. We conclude that A„(q:*) is a multiple eigenvalue. □ 

Appendix A. Discrete Graphs 

In this section we consider the analogues for discrete graphs of some of the theorems proved 


above, namely theorems 3.1, 3.3, and 3.4 


As for the mandarin graphs, one can consider their discrete analogues by placing several 
intermediate degree-two vertices per edge. However, the anomaly of the nodal count (theo¬ 


rem 


3.1) is only partially exhibited in this case. More specifically, in numerical experiments 
we saw that the nodal surplus stays anomalous (i.e., 7 ^ 0 and 7 ^ (3) only in the bottom 

half of the spectrum. Increasing the number of intermediate points per edge, one can ap¬ 
proximate any given eigenfunction of the quantum graph, but with a discrete eigenfunction 


that stays “low” in the spectrum, so these observations do not contradict theorem |3.1 


We introduce the relevant definitions for discrete graphs in section A.l and discuss the 
extrema of dispersion relations of inhnite periodic discrete graphs in section A.2 


(Hf)(u) = -Y,f(v) + q(u}f[u} 


A.l. Introduction to discrete graphs. Let L = {V,E) be a simple connected finite graph 
with vertex set V and edge set E. We define the Schrodinger operator with the potential 
g : 1 / -)■ M by 

(A.l) H : ^ 

That is, the matrix H is 

(A.2) H = Q-C, 

where Q is the diagonal matrix of site potentials q {u) and C is the adjacency matrix of the 
graph. It is perhaps more usual (and physically motivated) to represent the Hamiltonian 
as H = Q + L, where the Laplacian L is given hy L = D — C with D being the diagonal 
matrix of vertex degrees. Since we will not be imposing any restrictions on the potential Q, 
we absorb the matrix D into Q. The operator H has \V\ eigenvalues, which we number in 
increasing order as before. 

The nodal count of a (real) eigenfunction /„ is dehned as the number of edges on which 
the eigenfunction changes sign, i.e., 

{(M,n) e E : fn{u)fn{v) < 0} 

This count is most relevant for eigenfunctions which do not vanish at vertices. Otherwise, 
there exists alternative definitions |BLSn71 IBOSOSj . but they are not relevant for the results 
in this paper. 

We define the magnetic Hamiltonian (magnetic Schrodinger operator) on discrete graphs 
as 

(A.3) (H^f) (w) = - (^) + Q{u)f (u ), 


V'^U 
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Figure A.l. An infinite Z^-periodic graph with hexagonal lattice and chosen 
fundamental domain in box (a); the chosen fundamental domain with quasi- 
connected vertices marked with the same shape and color (b). 



with the convention that = —Au^v, which makes self-adjoint. For further details, 
the reader should consult |LL931 ISun941 ICdV981 ICdVTHTllj . 

A sequence of vertices C = [mi, U 2 , ■ ■ ■, Un] is called a cycle if each two consecutive vertices, 
Uj and Mj+i, are connected by an edge (wn+i is understood as Mi). The flux through the cycle 
C is defined as 


(A.4) 


— {^ui,u 2 + • • • + Au^_-^^u^ -j- Au^^ui) mod 27r. 


Two operators which have the same flux through every cycle C are unitarily equivalent (by 
a gauge transformation). Therefore, the effect of the magnetic held on the spectrum is fully 
determined by /3 huxes through a chosen set of basis cycles of the cycle space. We denote 
them hy ai,..., ay and consider the n-th eigenvalue of the graph as a function of ct. 

A result similar to theorem 2.3 holds for discrete graphs (in fact, discrete graphs were the 
original context in which the magnetic-nodal connection was established |Berl31 ICdV13] i. 
Note that in the discrete version of the theorem the nodal count should be modihed if one 
considers a symmetric point a* G {0,7r}^ which differs from zero. The nodal count one 
should consider in such a case is 




{{u,v) e E : {OL*) fn{u)fn{v) > 0} 


and its surplus, 0^ — (u — 1), equals the Morse index of An at the the symmetric point, ck*, 
according to the corresponding theorem |Berl31ICdV13j . Note that this modified nodal count 
is identical to the one previously defined if a* = 0. 


A.2. Periodic discrete graphs. Let G be an inhnite Z^-periodic discrete graph such that 
its Hamiltonian has a Z*^-periodic electric potential, q, and no magnetic potential. We denote 
by Sj{v) the vertex in G resulting from shifting the vertex v in the positive Z^-direction. 


We choose a fundamental domain W and dehne two vertices (u, u) G W to be the quasi- 
connected pair if the vertex Sj{v) is connected to u in G (see hgure A.l). In general, a pair 
of quasi-connected vertices is not connected in W. We restrict ourselves to graphs G that 
have only one quasi-connected pair in each direction, at least for some choice of W. 
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Let CK = {ai,a 2 ,- ■ ■, and on the fnndamental domain W define the operator = 
Q — C — M“, where Q is a diagonal matrix with real entries representing electric potential, 
C is the connectivity matrix of W, and 

if u'j is the quasi-connected pair, 
if y'j is the quasi-connected pair, 

0 otherwise. 


Mm = 


Similarly, let 7 = (do 72, 
M'*' is the diagonal matrix 


M? , = 

(v,v) 


., 7 fc) and on W define the operator iL'’' = Q — C — M'*' where 


if {v,u) is the quasi-connected pair, 
^ if (u, v) is the quasi-connected pair, 

o' 


otherwise. 

We are now ready to present the main theorem of this appendix, which is the discrete 


analogue of theorem 5.3 


Theorem A.l. Suppose the infinite periodic discrete graph G has no magnetic potential and 
only one guasi-connected vertex pair in each direction, A„(q:) has a critical point cx* that is 
not at a symmetry point of the Brillouin zone (i.e., such that a* 7 ^ 0 , 7 rj, the eigenvalue 
Xn{(x*) is simple, and the corresponding eigenvector f is non-zero at all guasi-connected 
vertex pairs. Then 

( 1 ) 7 * = is real for all j where {v,u) is the guasi-connected pair, 

(2) A = An(Q:*) is a degenerate eigenvalue of H'^*, and 

(3) if, additionally, W is a tree, there exists an internal vertex of W, of degree three or 
higher, such that the eigenfunction f is zero at this vertex. 


In order to prove theorem |A.l 


we establish the following two lemmas. 


Lemma A.2. Suppose that the discrete infinite periodic graph G has no magnetic potential 
and only one guasi-connected vertex pair in each direction. If cx* is a critical point 0 /A„(q:) 
and Xn{cx*) is a simple eigenvalue, then the eigenvector f corresponding to A„(q:*) satisfies 

(A.5) Vj = l,2,...,fc 

where {v,u) is the guasi-connected pair. 


Proof. Since An(cK*) is simple, Xn{(x) is analytic around the critical point cx = ck*. Let fa be 
an eigenvector of norm one corresponding to A„(q:). In particular, / = fa*. We have 

—ie^ocj jg jth quasi-connected pair, 

ie-'^^j if u) is the quasi-connected pair, 

0 otherwise. 


dH^ 
5a 1 


{v,u) 


dM^ 

da.i 


(VyU) 


Therefore, since cx* is a critical point of Xn{cx), one can see that 


d 

0 = ^An(Q:^ 


which completes the proof. 


5a,- 


/,/ = -fe*“^7(n)/(n) + ie *“^/(n)/(M) 


= 2 Im 




□ 
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Lemma A. 3. Suppose that G is a discrete infinite periodic graph and a* is a critical point 
of Xn{cx.). If the eigenvector f of H°'* corresponding to X = Xn{cx*) is non-zero at all guasi- 
connected vertex pairs, then f is also an eigenvector of H~* corresponding to the same eigen¬ 
value X where {v, u) is the guasi-connected vertex pair and 

7 ; = eM. 

f{v) 

Proof. We will demonstrate that H'^* f = H^* f = Xf. Using the definitions of the operators, 
one can see that this is equivalent to showing 

H^*f = {Q- C)f - M-^^f = {Q- C)f - M-*f = H^*f = Xf, 

or in other words 

(A.6) M^7 = M“7- 

Suppose that (u, u) is the quasi-connected vertex pair. Then at v we have 
(Ar7)(!.) = 7 -/(!') = = (M-'/Kv). 

Similarly, at u we have 

(M-'7)(«) = fm = e-‘-'.f(v) = (Af“7)(f). 

'Xj 

and at a vertex w which is not in a quasi-connected pair, 

□ 


Proof of theorem A.l. By lemma A.3, the eigenvector / of is also an eigenvector of H~^ 


(with the same eigenvalue A). Choose j such that a* 7 ^ 0 , tt (such a* exists by the theorem’s 


conditions). By lemma A.2 


e*“77)/(T) 


for all j, which means that / has some non-real entries and cannot be made real by scalar 
multiplication. However, it also means that the operator is real and self-adjoint. There¬ 
fore, the real and imaginary parts of / must both be linearly independent eigenvectors of 
H~^*, which implies that A is a degenerate eigenvalue of if'*'*. 

If W is a tree, a degenerate eigenvalue can only occur if there exists an internal vertex of 
degree at least three at which all eigenvectors from the eigenspace vanish |Fie75] . □ 
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